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ABSTRACT 


Physical-and-mathematical  model  of  interaction  processes  in  cumulative  jet  (CJ)  influenced  by 
pulsed  electric  current  as  physical  basis  of  electrodynamics  armor  (EDA)  has  been  studied. 
Magnetohydrodynamics  model  of  CJ  with  current  has  been  formulated.  The  role  of  the  current 
integral  and  the  rate  of  the  discharge  current  growth  as  the  main  factors  influencing  the  penetration 
deepness  decrease  have  been  shown.  The  critical  magnitudes  the  both  parameters  as  dependencies 
on  a  jet  charge  caliber  has  been  calculated.  The  parametric  analysis  of  coordination  of  the  discharge 
circuit  parameters  and  dimensions  of  EDA  on  the  base  of  numerical  modeling  has  been  done.  The 
theoretical  estimations  provide  the  design  of  EDA  experimental  assemblies  for  EDA  efficiency 
demonstration  and  the  theoretical  model  checking.  Two  groups  of  experiments  namely:  with 
laboratory  jet  charges  and  with  standard  80-mm  jet  charges  has  been  conducted.  The  experimental 
and  theoretical  results  are  discussed.  Technical  problems  of  EDA  realization  and  experimental 
modeling  are  discussed.  It  has  been  demonstrated  that  two  and  more  times’  decrease  of  the 
penetrating  deepness  can  be  provided  at  battery  capacity  of  units  of  mF,  discharge  voltage  5-10  kV 
for  jet  charge  caliber  20  -  80  mm. 
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ABBREVIATION  AND  DESIGNATION 


Explosive  substance 
Capacitor  battery 
Cumulative  jet 
Electrodynamics  armor 
Electromotive  force 
CB  energy 

Energy  allocated  in  active  resistance 
CB  voltage 
Electric  current 
Discharge  current 
Discharge  current  maximum  value 
Time  of  the  J 0  -  J m  achievement 
Discharge  current  average  value 
CB  charge 

Discharge  current  critical  value 
Critical  maximum  value  of  the  current  in 
the  cavity 

Time  of  J cr  current  formation 
Magnetic  field  induction 
CB  capacity 

Interelectrode  spacing  value 
Volume  element 
Active  part  of  CJ  mass 
caliber 

Electrical  field 

Average  density  of  the  current  in  CJ 
Density  of  induced  current 
Circuit  parameter 

Current  integral  for  conductor  explosive 
vaporization 

Current  integral  in  electrodynamics  model 
Atom  ionization  potential 

Jet  part  length  between  electrode  and 
contact  surface 

CJ  penetration  deepness  to  the  target 
CJ-without-current-penetration  deepness 
CJ-with-current-penetration  Deepness 
Circuit  inductance 
Circuit  inductance  per  unit  of  length 
Atom  concentration  in  plasma 
Atom  concentration  in  metal 
Pressure 
Critical  pressure 


PG  Hydrodynamics  pressure 

PM  Magnetic  pressure 

r  CJ  cross-section  radius 

r  CJ  resistance  average  value 

r0  CJ  initial  radius 

Tj  Plasma  channel  radius 

rc  Cavity  radius 

rs ,  6  Contact  Surface  Spherical  Coordinates 
y/  Material  strength  tensile 

R  Circuit  active  resistance 

R  CJ  external  radius 

R0  Resistance  value  per  unit  of  length 

Rc  CJ  resistance  per  unit  of  length 

Rp  Plasma  channel  resistance 

S  CJ  cross-section 

t  Penetration  time 

Tbh  Full  time  of  CJ  penetration 

Tcr  Critical  temperature 

TF  Fermi  temperature 

Tj  Current  pulse  duration 

Tv  Atom  vaporization  temperature 

U  CJ  penetration  rate  with  current 

U0  CJ  penetration  rate  without  current 

All  Penetration  rate  valuation 

V (Vr ,  V z )  layer  velocity,  C J  velocity 
V0  Initial  CJ  velocity 

V  Rod  critical  rate 

cr 

V p  Plasma  rate 

Y  Material  dynamic  hardness 

Z  Penetration  deepness  current  value 

a  Angle  between  CJ  and  normal  to 

electrode 

8  Damping  decrement 

8ck  Skin-layer  thickness 

p  Magnetic  permeability 

p  electrode  material  density 


pc  Jet  material  density 

cr  CJ  material  conductivity 

<70  Plasma  conductivity 

v  CJ  material  viscosity 


Plasma  velocity  to  penetration  rate  ratio 
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INTRODUCTION 


As  it  has  been  shown  theoretically  in  [1,  2]  cumulative  jet  (CJ)  penetrability  can  be  reduced  by 
passing  a  pulsed  current  through  CJ.  The  numerical  estimation  showed  that  this  principal  idea  might  be  used 
for  the  realization  of  simple  in  design  and  technology  electrodynamics  armor  (EDA)  against  jet  charges  (JC). 
The  technology  of  demonstration  experiments  also  could  be  easy  realized.  On  the  other  hand  the 
development  of  a  reliable  EDA  using  so  limited  means  as  the  discharge  of  a  capacitor  battery  (CB)  through 
CJ  does  necessary  to  provide  a  number  of  modeling  (methodical)  experiments.  This  contract  with  limited 
financing  did  not  give  possibilities  to  provide  “full-scale”  research  on  all  scientific  and  technological  aspects 
of  EDA.  The  main  task  of  the  present  study  is  to  develop  a  physical  model  of  processes  in  EDA  and  to 
demonstrate  the  CJ  penetrating  deepness  decrease  under  high  pulsed  current  influence. 

The  presented  report  adds  up  the  results  of  about  one  year  activity  of  R&D  PLATONIK  scientific 


group. 

In  the  report  the  parametric  analysis  of  experimental  results  on  cumulative  jet  (CJ)  effected  by  a 
pulsed  current  is  carried  out. 

The  parametric  analysis  of  theoretical  and  experimental  results  puts  the  following  problems: 

1.  To  show  that  the  offered  discharge  physical-and-mathematical  model  takes  into  account  basic  physical 
processes  in  CJ  under  effect  on  a  pulsed  current.  It  allows  to  predict  with  an  acceptable  error  relation  of 
CJ  penetration  deepness  in  EDA  from  main  recorded  parameters  -  capacity  and  voltage  of  the  condenser 
battery  (CB),  feeder  line  properties  connecting  CB  and  EDA  in  view  of  the  CJ  actual  parameters  (lengths 
and  radius,  maximum  velocity,  material).  The  analysis  should  takes  into  account  the  EDA  structural 
diagram  (quantity  of  electrodes,  distance  between  them)  and  protected  wall  material  performance. 

2.  To  select  among  the  current  pulse  parameters  (steepness  and  duration,  maximum  and  average  discharge 
current  values)  the  main  ones  determining  the  EDA  efficiency. 

3.  To  determine  influence  some  parameters  not  arrested  in  everyone  experiment,  for  example,  jet  resistance 

and  resistance  of  "jet  -  electrode"  contact  on  EDA  efficiency,  to  show  necessity  to  take  into  account 
them  at  statistical  processing  of  the  experiment  results,  to  formulate  requirements  to  decreasing  their 
influence  on  EDA  efficiency. 

The  physical-and-mathematical  model  of  electrodynamics  armor  (EDA)  that  allows  systematizing 
the  experimental  results  and  opens  a  capability  to  forecast  the  process  parameters  carries  out  such  analysis. 

The  present  report  consists  of  3  chapters  and  4  appendices. 

In  chapter  1  the  design  formulas  for  CJ  with  current  penetration  deepness  determination  in  different 
models  described  in  Appendices  1  ...  3  are  presented.  The  determination  scheme  of  average  current  through 
CJ  at  given  discharge  circuit  performance  is  shown. 

In  chapter  2  the  laboratory  experimental  researches  on  CJ  penetration  into  electrodynamics  bumper 
system  are  described.  The  penetration  deepness  is  compared  with  calculating  value  depending  on  circuit 


R 

parameter  y  =  — 


[C 

I —  ,  CB  capacity  and  energy. 

V  Li 
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In  chapter  3  the  scheme  and  results  of  EDA  full-scale  tests  by  usage  of  nominal  cumulative  means 
are  presented.  The  CJ  of  90  mm  caliber  penetration  computational  deepness  as  a  function  of  the  CB  capacity 
and  energy  for  different  circuit  resistance  are  presented.  Appendices  1  ...  4  represent  the  computational 
models  for  CJ  with  current  penetration  deepness  estimation. 

In  the  report  there  are  conclusions  by  major  experimental  and  computational  research  results. 

This  report  is  final  and  is  made  according  to  item  0003  of  Contract  No.  68171-98-M-5655 
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1  ELECTRODYNAMICS  ARMOR  PARAMETERS  COORDINATION 
_ _ PROBLEM _ 

EDA  parameter  coordination  problem  poses  as  a  problem  of  its  parametric  optimization  for  the 
purpose  of  EDA  maximum  efficiency.  EDA  efficiency  is  estimated  by  CJ  penetration  decreasing  factor 

>1  (L0  -  penetration  deepness  without  current,  Lm  -  penetration  deepness  with  current).  Values  L0,Lm 

Lm 

are  determined  from  the  cavity  axis  dimension  in  a  target  (monolith  bumper)  in  which  the  electrodes 
connected  with  CB. 

The  EDA  parameters  can  be  the  following: 

•  adjusting  (capacity  C,  voltage  Vc)\ 

•  defining  by  the  feeder  line  scheme  selection  and  electrode  system  (discharge  circuit  inductance); 

•  depending  on  current  closing  conditions  in  interelectrode  spacing  and  not  inspected  in  experiment  (active 
discharge  circuit  resistance  R ,  including  resistance  of  a  discharge  “electrode  -  jet”  gap). 

As  a  result  the  solution  of  the  EDA  parameter  coordination  problem  the  at  which  the  influence  of  the 
last  uncontrolled  factor  would  be  minimum  condition  should  be  determined.  The  problem  is  decided  by 
usage  of  the  physical-and-mathematical  model  permitting  to  predict  the  EDA  efficiency  function  from  its 
parameters. 

The  EDA  parameter  coordination  problem  breaks  up  on  some  individual  problems. 

1. 1  CONSEQUENCES  FROM  PHYSICAL  EDA  MODEL 


The  basis  of  a  physical  model  of  the  electrodynamics  interaction  on  CJ  is  the  description  of  the 


fjJ  0 

generation  mechanism  of  the  induced  by  CJ  motion  currents  jr  ~  oVcB  (B  —  —  —  magnetic  field 

270*0 

induction  created  by  discharge  current  J0 ;  V0,r0  -  CJ  velocity  and  radius,  cr  -  CJ  material  conductivity)  in 
discharge  current  magnetic  field.  In  case  of  pulse  current  interaction  the  Lorentz  force  braking  a  material  is 

i 

affixed  on  a  thin  CJ  skin-layer  ((  -  current  pulse  duration,  p  -  magnetic  constant,  8CK 

thermal  energy  generated  per  unit  of  volume  in  time  unit  is  ~  oB2Vj . 

Jet  kinetic  energy  converted  to  thermal  one  for  the  time  8t  is  calculated  as: 

oB2v02a  ~ 


j  TKSjJ. 


-).  The 


( p  -  CJ  material  density). 

Taking  into  account  cr  ~  106  (nxm)~‘  ,  p  =  8x!03  kg/m3  it  is  possible  to  receive  J  ~  100  L4, 
registered  in  the  electrodynamics  interaction  on  CJ  experiments  [1,  2]  at  St  ~  10  3  sec. 
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The  specific  heat  of  copper  vaporization  at  temperature  is  higher  than  melting  point  [3]  that  is  equals 

V2 

to  4. 7xl06  J/kg.  The  specific  jet  kinetic  energy  is  - 3xl07  J/kg.  It  is  possible  to  expect  an  explosive 

evaporation  of  CJ  mass  concluded  in  a  skin-layer  if  the  specific  thermal  energy  generated  in  a  skin-layer 
exceeds  ~  20  %  of  jet  specific  kinetic  energy.  In  this  case  the  plasma  sheath  can  be  organized  near  to  a  jet. 
Thus,  the  plasma  located  in  high-pressure  area  will  be  possible  to  have  imperfect  plasma  properties. 

The  latter  is  proved  by  estimation  of  the  pressure  and  temperature  critical  values  determining  plasma 
transition  in  the  condition  with  metallic  conduction.  For  copper  this  estimation  can  be  conducted  under  the 
formula: 

K  =0.411;  T.  =0.041  ion, 

where: 

/.  =  7.5  eV. 

ion 

Whence  Pcr  ~  10s  N/m2,  Tcr  ~  0.3  eV. 

Definite  thus  pressure  is  about  magnetic  field  pressure  (at  the  current  of  about  ~ 100  kA): 


P  =  -^-~ 
M  8K2r2 


10s  N/m2 


and  is  less  then  jet  hydrodynamics  pressure: 


/>  =  ZJL 


pVo 


10“  N/m2 


The  plasma  channel  which  conductivity  can  be  close  to  a  conductivity  of  a  metal  can  short  current.  It 
could  be  if  a  current  in  metallic  conductor  flows  inside  thick  skin-layer  which  thickness  is  8ck  «  r0  So, 
according  to  non-classical  skin-layer  model  for  the  plasma  channel  [1]  it  is  possible  to  consider  that  the 
magnetic  field  and  the  current  occupies  the  whole  plasma  volume. 

As  both  plasma  critical  temperature  Tcr  ~ 0.3  eV ~  4x1  03  K  and  temperature  relevant  to  latent  heat 

2 

of  atom  vaporization  ( Tv  ~  2  eV  ~  3x10*  K)  are  less  than  Fermi  temperature  ( TF  =  4.3  x  10~“n3  ~  10s  K 
at  n  ~  1023  cm3)  the  plasma  conductivity  can  be  determined  as  in  [3]: 


cr  ~  <J  n 


n 


nn 


where: 

o0  =  JO6  (fixm)'1; 

n,  n0  atom  concentration  in  plasma  and  metal,  correspondingly. 

1  y2 

As  earlier  [1,  2]  the  model  in  which  n - —  and  a  =  cr0  —r-  is  used. 
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The  linear  resistance  R  of  the  plasma  channel  at  r0  <  r  <  r,  corresponds  to  the  jet  cross- 

section: 

1  n  r 

—  =  |  cr{r)2nrdr  =  cr0  2^r/  In  — . 

Rp  r0  ro 

The  jet  linear  resistance  Rc  depends  on  the  jet  radius  and  the  skin-layer  thickness  SCK  as 

follows: 

— —  =  <7C  2nr 0SCK ;  SCK  ~ 

R-c 

Taking  into  account  a0  ~  106  (fixm)'1,  ac  *  4. 7  x  106  (nxm)~!  (copper  conductivity  at  melt 

■p 

point),  SCK  (at  r  ~  5  nsec)  ~  0.5 mm,  r0~  1mm,  In—  «  2  it  is  possible  to  receive 

r0 

Rp  ~Rc~0.35xl0-2  fi 

Thus,  the  estimations  give  the  identical  value  for  the  resistance  of  both  the  plasma  channel  and  the 
jet.  The  interelectrode  spacing  resistance  might  be  determined  experimentally  and  being  presented  in  Fig.  1  is 
close  to  the  obtained  estimation. 

In  a  plasma  flowing  from  the  cavity  with  the  velocity  about  Vp  an  EMF  can  be  generated  as  a  result 
of  the  interaction  between  discharge  current  J0  and  a  magnetic  field  in  the  cavity  [1]. 

EMF  =  \BV  dr  =  V  ^ In ^  =  V  J0Li0, 

J  p  p  2n  r, 

ri  1 

where: 

Li0  linear  inductance  of  a  current  circuit  in  the  cavity. 

When  the  plasma  flow  velocity  achieves  a  value  of  rate  V p  a  loop  with  a  current  that  is  irrelevant  to 

the  discharge  current  can  appear  in  the  cavity.  This  model  assumes  that  the  energy  of  the  induced  current  in 
the  cavity  is  taken  mainly  from  the  CJ  kinetic  energy  but  not  from  CB  storage  energy. 

The  implementation  of  the  following  scheme  of  CJ  penetration  process  is  possible  under  the  electric 
current  passing.  At  CJ  velocity  5  ...  10  Icm/s  in  a  result  of  CJ  partly  skin-layer  evaporation  in  the  cavity  a 
plasma  can  appear.  In  the  magnetic  field  of  the  discharge  current  in  the  plasma  the  current  of  density 
jr  ~  aVpB  is  induced.  The  interaction  of  the  induced  current  with  the  magnetic  field  the  Lorentz  forces 

braking  the  flow  arises.  The  pressure  generating  in  the  plasma  is  transmitted  to  CJ  through  the  contact 
surface  of  the  CJ  head  in  the  cavity.  This  process  results  decrease  of  the  penetration  velocity. 

The  mechanism  of  jet  kinetic  energy  losses  is  possible  at  the  radial  target  deformation  caused  by  the 
CJ  impact.  As  a  result  of  radial  velocity  Vr  the  induced  current  of  density  jz  ~  aV rB  appears.  Both 
discharge  and  induced  currents  has  to  be  similar  direction.  In  this  case  the  braking  force  operats  against  the 
magnetic  pressure  determined  by  current  J 0  flowing  in  CJ  or  in  the  adjacent  plasma  channel.  The  magnetic 
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pressure  at  J0  ~  100  kA  is  equal  to  PM  =  ~  108  N/m2.  It  is  less  than  the  hydrodynamics  pressure 

8n  r0 


component  that  equals  to  Pa 


10"  N/m2.  However,  the  magnetic  pressure  force  is  determined  as  a 


product  of  the  pressure  on  CJ  lateral  surface  area,  meanwhile  the  hydrodynamics  pressure  force  is 
determined  by  CJ  cross-section.  Therefore  at  axial  current  through  CJ  the  essential  losses  of  the 
hydrodynamics  pressure  are  possible.  It  offers  to  consider  two  mechanisms  of  electrodynamics  action  on  CJ. 
The  first  one  actuates  plasma  deceleration  by  the  magnetic  field.  In  so  doing,  pressure  pp{r)  in  the  plasma 


appears.  Then  CJ  penetration  rate  U  in  the  target  (contact  surface  movement  velocity)  can  be  estimated 
from  Bernoulli's  equation: 


P  = 


pu 2 


+  V  +  Pp  = 


Pc(y.-uf 


2 

J  pp(r)cos$dS 


dS  =  2  nr 2  sinOdO 


nrn 


(1) 


where: 

p ,  pc  target  and  CJ  material  densities; 
y/  target  metal  strength  tensile; 

V0  CJ  initial  velocity; 

dS  element  of  the  contact  surface  (written  in  spherical  coordinate  system); 

rs,9  spherical  coordinates. 

From  Bernoulli’s  equation  (1)  the  pressure  Pp  effects  on  CJ  similar  as  the  target  material  resistance 

defined  by  its  strength  tensile  y/ .  It  leads  to  CJ  penetration  velocity  U  decreases  (see  Appendix  4).  On  the 
other  hand  the  additional  pressure  onto  the  contact  surface  has  to  cause  the  cavity  diameter  increase  as  it  is 
watched  in  experiments. 

The  second  mechanism  is  caused  by  the  hydrodynamics  pressure  action  on  the  contact  surface  at  the 
radial  deformation  in  the  magnetic  field.  For  non-steady-state  flow  Bernoulli’s  equation  is  used  for 
qualitative  estimation  of  CJ  velocity  influence  on  the  penetration  rate.  If  to  neglect  the  target  material  elastic 

resistance  and  to  consider  p~  pc  it  is  possible  to  receive: 

pu!(t)  PcW)-u(i) r 

2  2 

whence: 

For  this  mechanism  CJ  penetration  velocity  decrease  is  explained  by  CJ  electrodynamics  braking  (jet 
velocity  decrease). 
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1.2  INDUCED  CURRENTS  IN  CUMULATIVE  JET 

In  [1,  2]  it  has  been  shown,  that  the  EDA  efficiency  depends  on  currents  induced  in  the  cavity.  The 
model  assumes  the  formation  of  current  loops  in  the  cavity  irrelevant  with  the  discharge  current  circuit.  Such 

dJ0 

currents  can  not  proceed  the  critical  value  defined  by  the  discharge  current  building  up  : 

J°r^~dttcr'  (2) 

where: 

tcr  time  of  the  current  growth  in  the  front  of  the  discharge  current  pulse. 

Fig.  1  represents  the  experimental  functions  of  discharge  current  J  0(t)  and  resistance  R(t)  of  t  at 
the  conditions:  CB  capacity  C~6xJ0  3  F,  CB  voltages  Vc  —  21kV,  CC  caliber  D  — 100  mm.  The 

interelectrode  spacing  was  executed  under  the  double-electrode  scheme. 

Fig.  1  represents  -  .  . 

l,MA  R.mn 

also  the  computational  1 4  _ ___ _ _ _ _ _ 

current  pulse  J  =  Jcr.  1,2-12  — - J0  (experiment) 

The  factor  of  the  1  10  yX  \ 

penetration  deepness  0,8  ■  8  /  \ 

/  R  ( experiment j  \  I 

decrease  for  the  given  06.  6  /  /\  /~\.  \ 

pulse  duration  Tj  can  be  0,4  -  4  -y/  . \ 

calculated  from  the  0,2  -  \ 

formula  A1 .5  0  - « - - - 1 - — 1  1  ' 

0  10  20  30  40  60  60  70  80 

(Appendix  1):  T.ps  I 

Fig.  1  Experimental  and  computational  pulsed  currents  through  CJ  in  the 

cavity  (test  001) 


•O  v - 

=S-  =I  +  ^-,A0  = 

LJf,  A° 


6nM0  _  dJ0  ~  9  j  -  r  T 

J—-,Jcr-7  -lcr’Y  *  ’  1 1  Jcr1i 

cr0JuV0  dt 


where: 


M0  =  pnr20l0 
cr0  =  106  (Qxm) 


current  integral; 

mass  of  CJ  active  part  with  length  of  l0 ; 


06  (Qxm)  1  plasma  conductivity; 

For  CJ  of  caliber  D  the  following  approximating  relations  were  taken: 


•  radius  r0 - D ; 
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•  total  penetration  deepness  L0  ~  4D ; 

•  jet  length  l0  ~  8D ; 

_  4 

•  total  penetration  time  without  current  IBH  ~  — . 

' o 

For  mode  presented  in  Fig.  1  TBH>Ti  (TBH  ~  100  psec,  Tt  ~  60  /usee)  during  period  of  time 
x  ~  TBtt  -  T  the  current  influence  may  be  neglect.  In  this  case  the  total  penetration  time  may  be  estimated 

BH  i 

by  the  formula  A1 .6  (Appendix  1): 


Fig.  2  gives  the  current  in  CJ  of  D  =  100  mm  for  CB  capacity  12.8  mF  (test  002)  and  19.2  mF  (test 


003),  at  voltages  5  kV. 

Initial  current  growth  rates  in  these  modes  are  different: 
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Lt  (003) 

z,  (002) 


2. 


For  CJ  of  100-mm  caliber  tcr  ~  10  jusec  (Appendix  2)  current  J0  <  Jcr{tcr )  in  tests  002  and  003 

is  observed  at  the  greater  part  of  the  discharge  pulses. 

In  this  case  the  current  estimation  in  cavity  at  average  value  of  the  CB  discharge  current  J  0  is 

possible  to  receive  as  (Appendix  2): 


( 


1-exp 


f  tW 


(5) 


Current  integral  equals  to: 


1 2  —  j ’  J 1  (t)dt  ~  yJ0tc 


\  1  c  J  J 


T 

l  +  ln— 

tcr  J 


(6) 


It  is  possible  to  accept  I  =  min{l  v 1 2)  as  an  estimated  value  of  the  current  integral. 


1.3  DETERMINATION  OF  DISCHARGE  CIRCUIT  PARAMETERS 

In  tests  001...  003  the  maximum  discharge  current  decreases  at  CB  capacity  increase  from  C  ~6 
mF  (test  001)  to  12  mF  (test  002).  In  test  003  (C  ~  18  mF)  maximum  discharge  current  stayed  on  the  level 
of  test  002.  Let  us  show  that  it  is  connected  with  the  change  of  the  active  resistance  of  the  discharge  circuit. 
Let  us  estimate  discharge  circuit  parameters  not  indicated  in  these  tests  using  CB  parameters  (C,Vc)  and 


current  curve  parameters  (maximum  current  value  J m  and  time  tm  of  its  achievement). 


giving  available  to  calculate  circuit  resistance  R  and  circuit  induction  Li , 
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T  =  2kJl£-, 


maximum  current  value; 


Tj  pulse  duration. 

Parameter  y  determination  allows  to  define  circuit  parameters  using  the  data  presented  in  Fig.  3  by 


calculating  the  following  criteria: 


y  = 


2  Jjm 

7t  CVc 


JJ 


m 


The  operational  sequences  for  y,LiyR  determination  are  presented  below. 
Test  001: 


y  =  0.26 ,  y  =  0.6,  =  0.45 ,  J00  =  Vc  y-  =  2.7xl06  A,  L  =  4xl0~7  Henry, 

J oo  V  A 


R  =  2yJ^  =  10-2  a 


The  calculated  value  R  is  close  to  measured  value  given  in  Fig.  1. 
Test  002: 


y  =  0.07 ,  ^-  =  0.3,  J 00  =  1.2  MA,  L  =  2.1xl0  7  Henry,  y  =  3,  R  =  24  mQ. 

Joo 

Test  003 : 

y  =  0.06  ,  —  =  0.3,  J00  = 1.2  MA,  L  =  4  x  J 0"7 Henry,  y  =  3,  R  =  27  mQ. 

J  00 

Thus,  the  essential  difference  of  discharge  parameters  in  described  tests  is  explained  by  parameter  y 
change  that  causes  the  change  of  the  circuit  resistance:  R(00l)~  10  mQ  R(002)~  R(003)~  25  mQ 


1.4  MODEL  OF  ELECTRODYNAMICS  EFFECT  ON  CUMULATIVE  JET 

The  theoretical  model  described  in  section  1.2  assumes  the  increase  of  target  resistance  to  CJ 
penetrating  in  the  case  of  plasma  generation  in  the  cavity.  It  is  supposed  that  the  plasma  is  generated  at  CJ 
skin-layer  explosion  evaporation  accompanied  by  the  CJ  electrodynamics  braking  in  the  magnetic  field  of 
the  discharge  current.  This  model  can  be  used  for  EDA  efficiency  estimation  when  the  discharge  current 
exceeds  a  critical  value  for  a  given  CJ  caliber.  At  smaller  currents  it  is  possible  to  guess  that  the  plasma 
influence  on  EDA  efficiency  is  negligible.  For  this  case  the  effect  of  the  penetration  deepness  decrease  might 
be  connected  with  CJ  kinetic  energy  losses  under  CJ  radial  deformation  in  the  magnetic  field  of  the  current 
passing  through  CJ.  CJ  deepness  penetration  can  be  estimated  from  the  analytical  solution  of  Navier-Stokes 
equations  (Appendix  3,  4).  At  different  assumptions  two  models  are  obtained. 


MODEL  I:. 
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Penetration  decrease  factor  is  obtained  as: 


2 

ft? 

1 _ 

1 - 

1 _ 

1 _ 

where: 


PM  = 


B2  juJ2 


2p  8jt2rl 

2 


P  -pKo 

rG  -  „ 


,(for  L0/Lm>l,  see  App.4) 


CJ  penetration  deepness  with  and  without  current; 
magnetic  pressure; 

hydrodynamic  pressure. 


(7) 


MODEL  2. 

Penetration  decrease  factor  is  obtained  as: 


-'0  _ 


1—  “1 

2 

r~  - | 

4 

3 

1 

1 

_ 1 

( for  Lo/Lm>l,  see  App.3  ) 


(8) 


In  these  models  the  EDA  efficiency  is  determined  by  two  parameters:  —  -  relation  between 

r0 

p 

penetration  deepness  and  CD  radius  when  Jo= 0,  and  — —  -  relation  between  magnetic  and  hydrodynamic 

Bg 

pressure  components. 

In  sections  3,  4  EDA  efficiency  estimations  on  different  models  and  the  comparison  with 
experimental  results  will  be  conducted. 

1.5  DETERMINATION  OF  DISCHARGE  CURRENT  AVERAGE  VALUE 

In  offered  models  EDA  efficiency  was  estimated  taking  the  average  (on  pulse  duration)  current 
value.  As  it  is  shown  in  section  1 .2  the  mode  with  maximum  efficiency  at  minimum  (critical)  capacity  Ccr 

can  be  determined  assuming  y  ~  y cr 


From  (9)  Ccr  is  defined  by  circuit  parameters  R}Lt . 
The  current  at  y  =  ycr  is  calculated  by  formula: 

f 


(9) 


Vc  (  R 


Jo  ~-j-texp 


\  2L>  J 


(10) 


The  average  current  value  in  the  critical  mode  is  (T;  >  ^bh  )  ■ 

j  Tbh 

J" =■=-  K(4* 


lBH  0 


V  c 

y  C^cr 

T 

J.  D  U 


(ID 
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The  average  current  value  in  the  mode  with  y  *  1  is  determined  by  formula  A2.7  (Appendix  2): 

.2 


j=j:  1+r  Tbk- 


V  Tt 


(12) 


-  2CV 

From  (12)  at  y  ~  0.8  the  average  current  value  is  J  = — —  This  value  corresponds  to  the 


Tj 


average  discharge  current  value  when  the  attenuation  can  be  neglected: 

-_q0  CVc  2CVc 


J  = 


T/  T/  T. 
/4  /4 


Equations  (10)-(12)  are  received  for  two-electrode-EDA.  In  the  case  of  tree-electrode-EDA  the 
values  calculated  by  these  equations  must  be  divided  to  2. 

1.6  EDA  OF  “ABSOLUTE  EFFICIENCY” 

For  estimation  it  can  be  determined  the  parameters  of  ideal  EDA  (for  example  assuming  factor 

^  ~  100 ). 


According  to  (A3. 13) 


^  =  2 
L. 


I 


P,  ’ 


(13) 


where:  l  =  d  +  Lm,  d  -  interelectrode  spacing. 
In  case  of  /  =  dcr  (Lm  «  d)  one  can  receive 


w  =,  4l  \Us. 

der  r0J  r  i  ~  p 

*  rM 


(14) 


Value  of  dcr  determines  the  critical  (minimum)  spacing  between  electrodes  for  ideal  EDA 

p 

efficiency.  At  J  =  400  kA,  r0  ~  1mm,  -2- =  125,  interelectrode  spacing  dcr  >  100  mm.  This  estimation 


is  right  at  very  high  current  growth  rate  when  the  penetration  effect  can  be  neglected. 
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2  LABORATORY  EXPERIMENTAL  INVESTIGATION  OF  CUMULATIVE 

JET  PENETRATION 

Principal  diagram  of  the  laboratory  experiment  is  presented  in  Fig.  4. 

In  the  first  experimental  series  cumulative  charge  mock-ups  ( 1 )  with  copper  conical  funnels  of  25- 
mm  caliber  were  used.  Two-section  CB  of 
total  capacity  600  /jF  was  connected  with 
electrodes  ( 4 ,  6)  by  a  low-inductive  high- 
voltage  coaxial  cable  (10). 

Basic  plate  9  assembled  from 
homogeneous  aluminum  alloy  plates  had 
total  thickness  of  about  200  mm.  The 
thickness  of  fiberglass  insulated  plate  (5)  has 
the  thickness  Aj  in  the  range  of  0...15  mm, 
and  the  thickness  of  plate  7  A2=5 ...  10  mm. 

CJ  penetration  deepness  without  current  was 
up  to  L0 uax  —  95  mm. 

Performance  of  three  tests  is 
presented  in  Table  1 . 

Current  oscillograms  registered  using 
/?C-integrator  output  connected  with 
inductive  sensor  11  are  given  in  Fig.  5. 


Table  1 


Test 

Calculation 

C,  mF 

V,  kV 

W,  kJ 

EXH 

L0,  mm 

Y 

warn* 

Lm,  mm 

01 

0.6 

■ 

5.6 

0 

25 

1 

13 

31 

3 

0? 

0.3 

2.8 

34 

0.2 

9 

0.15 

45 

2.1 

03 

wm 

60 

23 

0.25 

11 

34 

2.7 

Experimental  and  calculated  results  are  compared  in  table  1.  One  can  see  satisfactorily  agreement  of 
theoretically  predicted  data  and  experimental  ones.  The  data  were  calculated  at  J2  current  (Fig.  5).  For  test 

01  the  current  was  assumed  as  J  —  J  max - for  tests  02  and  03  the  average  value  of  the  discharge  current 

^  max 

was  used.  For  test  03  the  magnetic  pressure  was  taken  proportional  to  the  current  in  square  (calculated 
assuming  model  (Appendix  3)). 


C 


Fig.  4  Apparatus  scheme  for  EDA  laboratory 
investigation 

1-  Cumulative  charge,  2  -  Electric  detonator, 

3  -  Fiberglass  plate,  4,  6  -  Copper  electrodes  (1  mm),  5  - 
Fiberglass  (At,  mm),  7  -  Fiberglass  (A2,  mm),  8  -  Air 
Backlash  (3  mm),  9  -  Basic  plate  (target),  10  -  High- 
voltage  cable,  1 1  -  Inductive  sensor 
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The  formation  of  "ledge"  at  the 
cavity  inlet  for  case  of  a=60°  share  is  shown 
in  Fig.  6.  This  effect  apparently  is  connected 
with  the  action  of  electromagnetic  forces. 
However,  influence  of  angle  a  on 
penetration  deepness  as  it  follows  from 
calculation  results,  (table  1)  is  negligible. 

The  computational  estimation  of  CJ 
radial  expansion  in  the  cavity  for  test  01  is 
given  in  Fig.  7a.  Maximum  relative  CJ 


radius  increase 


is 


\roJ 


3.  It 


corresponds  to  7)  ~  Tm  condition.  CJ  and 
cavity  geometry  are  shown  in  Fig.  7  (not  in 
scale).  In  Fig.  7b  the  calculated  cavity  radius 
at  two-times  current  pulse  decrease  is  given. 
One  can  see  that  at  so  different  discharge 
parameters  essential  cavity  form  change  can 


6 

Fig.  6  Cavity  forms  for  a  25-mm 
cumulative  charge  in  aluminum 
target 

a)  Test  with  current  (03);  6) 
control  test  without  current 


J,  kA 


Test  01 


J!  -  current  in  upper  circuit  (Fig.  3),  J2  -  current  in  lower 
circuit 

Fig.  5  Current  function  of  time  in  tests  No.  10, 18, 19 


be  observed. 

Radial  deformation  for  the  case  of  CJ  without  current  is  seen 
only  in  place  round  the  contact  surface.  For  CJ  with  current  the 
deformation  is  spread  around  the  whole  cavity. 
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The  second  series  of  laboratory  experiments  was  carried  out  using  the  circuit  presented  in  Fig.  4,  and 
cumulative  charge  mock-ups  of  diameter  30,  50,  and  100  mm. 

Other  EDA  assemblies  parameters:  fiberglass  insulator  plates  (<5),  (7)  thickness  Ai=15  mm,  A2=10 
mm  (for  the  1st  experiment),  Ai=20  mm,  A2=15  mm  (for  the  2nd  experiment),  A;  =30  mm,  A2=20  mm  (for  the 
3rd  experiment).  CB  voltage  was  Vc  =  5  kV. 

CC  shooting  in  all  cases  was  normally  to  the  target. 

Experimental  CJ  penetration  deepness  L  and  energy  separated  in  active  circuit  resistance  Wr  as  a  part 
of  CB  energy  W  for  different  CJ  caliber  are  presented  in  Fig.  8. 

Calculated  parameters  for  100-mm  CJ  is  shown  in  Fig.  8. 

Power  losses  in  the  active  resistance  r(t)  for  r  is  equal  to 

wr(w)=]r(ty(t)dt. 

o 
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L,  mm  Wr,  kJ 
500  - 


K100  -  calculation 


Wr  > 

" 

K50 

j  CJ  current  can  be 

-  expressed  by  its  density: 

•/W  =  J  J(‘)dS=Sj(l). 

_  s 

'  Hence  it  is  possible  to  receive 

the  following  relation  for 

current  integral: 

r_  K(w) 

—  /=-^’  <U> 

where: 

r  average  in  CJ  cross- 

—  

120  section  S  value  of  the 
Wr,  Kflw 

resistance  at  moment  of  time 


Fig.  8  Relation  of  CJ  penetration  deepness  (L)  and  power 

dissipation  in  CJ  (Wr)  as  functions  of  CB  energy  (W).  K30,  K50,  It  ig  interesting  t0 

K100  are  charge  caliber.  The  experimental  results  of  IGSOAN  ^  current  integra| 

scientific  group  are  used _ J  ^  corresponding  to 

explosion  evaporation  (blasting)  of  a  conductor  with  observed  in  CJ  experiments.  For  example,  for  50-mm 
CC  at  W=40  kJ  and  Wr(W)=1.37  kJ,  CJ  diameter  ~2  mm  and  resistance  f  =  10~2  Q  the  current  integral  for  a 
conductor  blasting  is  estimated  as  ~1.95xlOn  A2xsecxm~4  [3]- 

In  the  described  experiments  average  value  of  the  current  integral  was  1=1.2  xlO1 6  A2 xsecxm4 . 

As  it  discussed  in  [1]  blasting  condition  in  the  CJ  skin-layer  has  not  to  depend  on  current  integral. 
This  condition  can  be  realized  at  the  achievement  of  some  critical  current  value.  For  example,  if 

r0  ~  1  mm  J cr  ~  200  kA. 

j  *  |  It  is  necessary  to  notice  that 

^  _ ^  ...  .1  in  the  model  of 

. *  — — “i  electrodynamics  interaction 

_ _ _  |  the  current  integral 

m 

T, 

a  ^  /  =  J  j(t)dt  determines  the 


a  b 

Fig.  9  Cavity  contours  for  a  50-mm  CC  and  aluminum  target, 
a)  Test  with  current,  b)  Control  test  without  current 


electric  charge,  passing 
through  CJ  for  the  pulse 


duration  Tt .  Current  integral  for  a  conductor  blasting  /  =  ^J2(t)dt  is  proportional  to  the  heating 

o 

power  density  exuding  in  conductor  body. 
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Fig.  9  illustrates  EDA  efficiency  in  comparison  of  two  cavity  pictures  drawn  in  the  same  scale  for  a 
50-mm  CC  with  current  (a)  and  without  current  (b)  in  aluminum  alloy  D16  target.  In  case  a)  at  current  J=430 
kA  L c=93  mm.  In  case  b)  L0=360  mm. 
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K100  equals  to  2. 

Fig.  10  represents  the  estimation  of  the  penetration  deepness  as  a  function  of  the  circuit  parameter 

y  =  —  /  /—  calculated  by  different  models  (Appendixes  1,  3,  4)  for  CJ  caliber  D  -  50  mm,  circuit 

2 

resistance  R  =  10~2  J Q,  and  inductance  Li  =  2  x  10~7  Henry. 

At  achievement  of  some  critical  value  y  ~  ycr  the  penetration  deepness  practically  does  not  depend 

on  y . 

The  calculated  penetration  deepness  as  a  function  of  the  CB  capacity  and  power  (at  CB  voltages 
Vc  =5x  103  V)  in  different  models  (Appendixes  1,  3,  4)  for  CJ  caliber  D  =  50  mm  is  presented  in  Fig.  1 1. 

CJ  penetration  deepness  as  a  function  of  CB  capacity  and  energy  is  given  in  Fig.  12  at  circuit 
inductance  Lt-4  xl  O'7  Henry. 

The  minimum  penetration  deepness  (Fig.  10,  Fig.  11)  does  not  depend  on  CB  capacity  (energy). 

CJ  penetration  deepness  as  a  functions  of  CB  capacity  and  energy  (at  voltages  Vc  =  5  x  JO3  V, 

R  ~  10~2  a  Li  =  2  x  10~7  Henry)  for  CC  of  different  caliber  (model  JVs  2,  Appendix  3)  are  given  in 
Fig.  13. 
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The  computational  and  experimental  (Fig.  8)  results  are  good  correlated  in  the  investigated  range  of 
parameters. 
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3  RESULTS  OF  EDA  FULL-SCALE  TESTS 


The  program  of  full-scale  tests  provides  the  demonstration  of  the  armor  efficiency  of  EDA  with 
capacitor-type  power  supply  when  full-scale  standard  cumulative  charges  are  used. 

The  tests  was  carried  out  by  a  scheme  presented  in  Fig.  14.  Bench  units  were  pointed  in  three  places. 
Steel  container  111  and  EDA  mock-up  2  are  placed  on  testing  pad  11  in  50  m  away  from  shelter  1.  CC  1  is 
established  near  the  shelter.  High-voltage  supply  unit  4  adjusted  from  control  panel  8  serves  for  CB  charging 
to  the  given  voltages  Vc .  CJ  penetration  initiates  CB  discharge.  Two  oscilloscopes  7  indicate  the  current 
pulse  by  registration  signals  of  inductive  sensors  9. 


Fig.  14.  EDA  testing  bench  diagram 

I  -  Shelter,  II  -  Testing  pad.  III  -  Steel  container,  1  -  Cumulative  charge,  2  -  Electrical  board, 
3  -  Capacitor  battery,  4  -  High-voltage  supply  unit,  5  -  Discharge  resistance,  6  -  Switch, 

7  -  Oscilloscope,  8  -  Control  panel,  9  -  Inductive  sensors,  10  -  Switch,  1 1  -  EDA 


Three-electrode  EDA  design  diagram  is  given  in  Fig.  15.  Outside  plate  1  is  added  mainly  for  the 
EDA  mock-up  protection  against  small-dimensional  striking  fragments.  The  basic  plate  was  executed  as  a 
package  of  armor  plates  of  total  thickness  250  mm. 

The  90  mm  cumulative  garnet  was  used.  The  range  of  CB  parameters  varying  was  the  following: 

•  l.OkJ  (C=1.0mF,  Vc=1.5  kV); 

•  5.0  kJ  (C=0. 4  mF,  Vc=5. 0  kV); 

•  12. 5  kJ  (C=l.  0  mF,  Vc=5. 0  kV) ; 

•  17.5  kJ  (C-1.4  mF  Vc=5. 0  kV). 
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Current  maximum  value  in  the 
experiments  achieved  300  kA  at  pulse  duration  20 
ju s,  the  minimum  current  was  80  kA  at  12  jus. 

CB  storage  energy  varied  in  1-17.5  kJ 
range  due  to  CB  capacity  varying. 

Average  penetration  deepness  Lo  in  case 
of  zero  current  was  determined  on  the  results  of 
measurements  in  five  independent  experiments.  It 
was  230  mm  at  deviation  within  190  mm...  250 
mm. 

Experimental  data  are  presented  in  Table  2 
and  in  Fig.  16. 


Fig.15.  Three-layer  EDA  configuration 


1  -  Armor  plate  (16  mm*),  2  -  Insulator  (30  mm), 

3,  6  -  Air  Clearance  (5  mm),  4,  7  -  Aluminum  plates 
(5  mm),  5  -  Insulator  (25  mm),  8  -  Basic  armor 
package  (target)  (220  mm) _ _ 

—  SSSSmSSSSSSSSSSSSSB -  "I  .ill  B  .i  ■■■»  gsOBBSS 

Table  2 


W,  kJ 

L 

AL 

LjL 

A >  /Ana* 

A  l^min 

1.0 

178 

48 

1.3 

1.5 

15 

5.0 

37 

59 

6.2 

46.0 

3.6  | 

12.5 

43 

76 

5.3 

57.0 

3.9 

17.5 

7.5 

7.0 

32.0 

48.0 

18.4 

In  Fig.  16  the  dependence  of  penetration  deepness  L  on  CB  energy  W  is  presented,  as  a  field  of 
experimental  points  not  united  by  a  curve. 


One  can  see  wide  spread  of  the 
results.  It  can  be  explained  by 
influence  a  number  of  factors. 

•  EDA  electrode  geometrical 
parameters  varying; 

•  lower  accuracy  of  the  penetration 
deepness  measurements; 

•  possible  dependence  of  the  circuit 
discharge  resistance  on  weather 
conditions. 

Never  less  one  can  see  that  at 
CB  energy  greater  than  5  kJ  the  EDA 
provides  more  than  five-time- 
decrease  of  the  penetration  deepness. 
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Fig.17.  Photos  of  the  Cavity  Braked  by 
Cumulative  Jet  in  Basic  Plate  (a  -  Without 
Current;  b  -  With  Current)  and  Input  Holes 
in  the  Bumper  System 


The  photos  of  the  cavity  in  basic  plate  (target) 
are  shown  in  Fig.17.  They  demonstrate  the  influence  of 
current  on  the  form  of  cavity.  The  cavity  in  the  absence 
of  the  current  (left  picture)  is  more  narrow  and  deeper. 
It  shows  that  under  the  current  enforce  the  dissipation 
of  CJ  kinetic  energy  take  place  in  a  larger  volume. 

The  computational  dependencies  of  the  CJ 
penetration  deepness  for  the  test  conditions  taking  in 
calculation  different  circuit  resistance  are  given  in  Fig. 
1 8  at  CC  caliber  90  mm. 

Circuit  discharge  resistance  can  be  variate  value 
mainly  in  the  force  of  the  equivocation  of  «jet  - 
electrode))  contact  resistance  that  practically  determines 
the  total  circuit  discharge  resistance.  The  results  of 
calculations  presented  in  Fig.  18,  19,  20  are  usable  for 
explanation  of  the  wide  spread  of  the  experimental  CJ 
penetration  deepness  data  shown  in  Fig.  16.  According 
to  Fig.  19  for  example  at  CB  energy  W=25  kJ  the 
resistance  valuation  in  two  times  causes  to  calculated 
penetration  of  the  detonator  operation  moment  deepness 
valuation  in  two  times. 

Other  reasons  for  the  experimental  data  spread 
could  be  the  dispersion  of  the  detonator  operation 
moment.  CJ  efficiency  decreasing  in  this  case  maybe 
estimated  as  use  of  CC  of  smaller  caliber. 
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COCNCLUSIONS 

1.  EDA  laboratory  static  tests  by  cumulative  means  of  25...  100  mm  caliber  shows  the  EDA  high  efficiency. 
Penetration  deepness  increasing  in  five  times  for  50  mm  CJ  achieved.  For  100  mm  caliber  this  value  is 
equals  to  three. 

2.  Dynamics  ground  test  result  with  90  mm  standard  cumulative  garnet  shows  high  instability  of  the  results. 
Maximum  decreasing  of  the  CJ  penetration  deepness  at  CB  power  of  about  25  kJ  is  higher  than  10  times. 

3.  At  laboratory  investigation  of  the  EDA  with  the  solid  rod  the  principal  possibility  of  its  efficiency 
decreasing  is  shown,  but  the  systematic  results  does  not  receive.  Penetration  deepness  decreasing  by  40  % 
of  the  7.62  mm  fairing  pin  was  achieved. 

4.  EMG  implementation  does  not  allow  to  receive  the  current  higher  than  350  kA,  because  of  the  generation 
and  EDA  corresponding  problem  did  not  solved. 

5.  The  physical-and-mathematical  process  in  EDA  model  usage  allows  to  systematized  the  experimental 
investigation  data  on  CJ  with  current  penetration  into  the  bumper  system.  The  possibility  of  its 
forecasting  in  EDA  in  wide  range  as  well  as  CJ  parameters  are  established. 

6.  The  minimum  CB  capacity  value  for  the  given  caliber  at  which  the  EDA  maximum  efficiency  achieved 
was  received  (CJ  minimum  penetration  deepness).  For  100  mm  CJ  at  CB  voltage  of  about  ~5  kV  the 
critical  CB  power  value  is  equals  to  20  kJ. 

7.  The  physical-and-mathematical  model  allowing  to  estimate  the  critical  (maximum)  value  of  the  discharge 
circuit  (interelectrode  spacing  resistance)  at  which  the  EDA  maximum  efficiency  achieved  in  wide  range 
of  CB  capacity  was  elaborated.  For  CJ  caliber  of  about  100  mm  this  resistance  is  equals  to  ~10  7  Q. 

8.  The  estimation  of  the  limited  (critical)  current  in  the  cavity  was  carried  out.  It  determines  by  the  leading 
edge  steepness  of  the  discharge  current  pulse. 
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APPENDIX  1 


ESTIMATION  OF  CRITICAL  CURRENT  AND  PENETRATION  DEEPNESS  VALUES 

A  CJ  motion  in  a  magnetic  field  causes  vortex  currents  in  the  CJ  skin-layer  providing  the  energy 
dissipation.  When  the  kinetic  energy  losses  become  comparable  with  CJ  bulk  material  evaporation  energy, 
the  explosive  evaporation  of  the  layer  is  to  be.  According  with  the  proposed  model  a  plasma  cloud  with 
metallic  conductivity  is  implemented  [1].  In  case  of  a  vortex-free  plasma  motion  (with  velocity  close  to  CJ 
speed)  in  the  magnetic  field  of  the  discharge  current  braking  Lorentz  forces  begin  to  act.  These  forces  are 
equilibrated  by  plasma  pressure.  The  pressure  in  the  plasma  is  transmitted  through  the  contact  surface  to  the 
CJ  head  resulting  the  penetration  rate  decrease. 

The  equation  of  an  energy  balance  in  the  cavity  [1]  is  the  following: 

'j]oV;B!U{t)2nrdrclt  =  -pnr20\vo -U(t))! -(V0  -£/,)J]  Au 

0rB  Z 

Taking  into  account  U0  -U(t)=  AU  and  expanding  (Al.l)  right  part  into  a  series  it  is 


possible  to  receive  the  approximate  relation: 

'fiaVf^p-rdrU(‘Mt  =  jp^-(K-UJsAU 


A1.2 


The  value  V0  -  U0  determines  the  jet  velocity  in  case  of  penetration  without  the  current  (in 
coordinate  system  in  which  the  contact  surface  is  fixed). 

In  this  case  the  penetration  rate  U0  defines  by  Bernoulli’s  integral  through  pressure  P  in 
braking  point  onto  contact  surface. 

rP-pul+Y=-2P0ty0-uJ,  ai.3 

where: 

Y  target  material  strength  tensile; 
p,  p0  target  and  CJ  materials  densities. 

Considering  that 

v,sr,-u,  ai.4 

from  (A1.2)  one  can  receive  (model  3): 


Where: 

L0,Lm  penetration  deepness  without  and  with  current  respectively; 


A1.5 


/  =  JcrTBH  current  integral; 
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CJ  penetration  time; 


_  6tM0 

*0  ~  2rr 

aju  V0 


M0  =  pnTg  1 0  CJ  active  part  mass. 

The  relation  (A1.5)  is  justice  when  penetration  time  TBH  is  less  than  current  pulse  time  7}  . 
In  case  TI  <  Tm 


1  ,  1 

- — +  7 - 

i  y 
y+T 


where: 

T 

y  —  —SIL  >  j 
Tj 

In  [1]  the  model  of  current/  excitation  in  the  cavity  by  EMF  deriving  in  the  plasma  layer  under  the 
influence  of  the  magnetic  field  on  moveable  plasma  flow  exhausting  out  of  the  cavity  has  been  considered. 
The  equation  for  current  J  calculation  has  the  following  form: 


R0ZJ  +  Li0Z  ^  +  Li0UJ  =  |  VpB(r)dr  =  VpJ0Li0 , 


where: 


_VJ o  . 


Z  penetration  deepness  current  value; 

J 0  discharge  current; 

Vp ,  U  plasma  velocity  and  penetration  rate; 

R0 ,  Lj0  CJ  linear  resistance  and  inductance.  _ 

Approximate  solution  of  (A  1.7)  in  case  when  the  discharge  current  determines  by  the  front 
edge  of  a  current  pulse  as: 


-t,Z  =  Ut, 


has  the  fo 


J=V^  dJo  Lw  1  |  exp{-y)  1  t 


U  dt  Rn 


^  U, 


As  - 2.,  V  =V0-U0,  y 
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Uy«,,J*Ly^Lt 


Al.  9 


It  is  assumed  for  plasma  conductivity  calculation: 

2 


<y  =  a, 


V 


A1.10 


where: 


r0  CJ  radius; 

r  varient  radius. 

Note:  According  to  [l ]  the  conductivity  of  the  imperfect  plasma  might  be  presented  in  the 

following  form  er  =  cr0  — ,  cr0  =  1 06  (jQ  x  m)~' ,  n0  ~  102s  sm~3 . 
no 

Plasma  resistance  R0  can  be  calculated  in  approximation  of  a  cable  line  including  external 
(rt  <r  <rc)  and  internal  (r0  <r  <r,)  conductors,  where  rc  -  cavity  radius. 

For  this  case 

j  rc  rc  r?  r  Al.  1 1 

—  =  f  cr{r)27trdr  =  cr0  \-\2vrdr  -  cr0  2 tit l  In— , 

Ro  J  i  r  r, 


tcr=^f-  =  M<r0 


Assuming  In—  ~  2  one  can  receive  the  critical  time  for  a  CJ  of  D  — 100  mm  a  value 


ter  ~  10  **■ 
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APPENDIX  2 


FORMATION  OF  A  CURRENT  CIRCUIT  IN  THE  CAVITY 

There  is  the  mode  when  current  J (see  AI.8  at  t  >  tcr)  does  not  achieve  its  critical  value  J  cr  in  the 
front  of  the  current  pulse.  In  this  case  the  current  in  the  closed  circuit  can  be  calculated  according  with  the 
following  assumption  J 0  =  J 0  =  const . 

In  this  case  the  solution  of  equation  A  1.1  has  the  form: 


J  =  Wc‘f 


f 


1-exp 


\  1  or  J  J 

‘cr 

I,  =  jj(t)dt  =  yJ0t 


,r=^~2 

'  u  . 


0  cr 


1 7  =  jj(t)dt  =  yJotcrlnj 


T, 


The  current  integral  is 


i  —  1  j  +  1 2  —  YJoh 


1  +  ln 


Tj 


*  cr  J 


A2.1 


A2.1 


Thus,  the  problem  of  the  determination  of  the  current  integral  reduces  to  the  problem  of  the 
calculation  of  the  average  discharge  current  J 0  . 

Analysis  of  the  experimental  data  shows  that  the  modes  with  CB  minimum  energy  when  the 
penetration  decrease  is  visible  are  close  to  the  critical  discharge  mode  {y  ~  1). 

In  this  case  Tj  ~  Tm  .  If  y  =  1  the  current  is  determined  as: 


j0=ftexp(-a). 


where: 


6  = 


R_ 

2L 


dumping  decrement. 


Let  us  determine  the  critical  value  of  the  discharge  current: 

T, 


_  y  /  m 

J?=T-T~  /'«?<-<»>*  = 


/  T 

-S  1 BH  0 


V  c 

-^,y  =  i- 


A2.3 


A2.2 


For  example,  if  V0  ~  5 x  103  V,  Li  -2x10  7 Henry,  R=  Iff2  Q,  Tj  =  Tm  ~  100 ps,  current  pulse 
duration  Tj(y  =  l)—>oo  (see  Fig.  3),  one  can  receive  J 0  =  400  kA. 
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To  estimate  J  0  we  need  to  define  the  corrected  relation  at  y  as  variate  (for  example,  at  CB 

capacity  valuation). 

It  gives  [3]: 

1.  in  case  y  <1 : 


coL 

2.  in  case  y  >  1 : 


J  =  —^-exp(-  St)Sina)t ,  co  -  co0 ^ J-y 2  , 


1  x  R  1  j?  C 

6>o=-7t=^’S  =— ’r  =  ~RJ— 


2L{  2  \L, 


A2.3 


o)*L 


A2A 


Expanding  in  series  Sin  cot  and  (e"*f  and  e~a*‘)  one  can  receive  the  approximate  relation  for 

current: 

1.  in  case  y  <  1 : 


J  =  -texp(-  Stil- 


L 


W 


A2.5 


2.  in  case  y  >  1 : 


J  =  -texp{-  <5n  1  + 


{co  *  t)2 


A2.6 


Averaging  the  relations  for  a  current  pulse  we  receive: 


7  _7 cr-l  +  y' 
J  0  —  J  0  *  2 


2r 4 

7  _7 crl  +  y2  Tl 


A2.7 


For  7]  <  Tbh  J 0  —  J 0  ^  2 


BH 


2y  T, 


At  7  =  7  J0=Jo- 

Relation  (A.2.7)  can  be  transformed  to: 


7  _VcCU  +  y2 
J  o  ~ 


T,  2  y' 


A2.8 


For  7  «  1  the  current  is  calculated  by  the  following  formula: 


where: 

q0  CB  capacity. 

Relation  (A2.8)  transmitted  to  relation  (A2.9)  at  y  ~0.8 . 


A2.9 
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APPENDIX  3 


MODEL  OF  THE  ELECTRODYNAMIC  BRAKING  UNDER  CJ  RADIAL 

DEFORMATION 


The  model  assumes  that  the  discharge  current  mainly  flows  through  CJ  and  closed  in  the  cavity  on 
the  contact  surface  separated  CJ  and  the  target.  The  effect  of  CJ  penetration  deepness  decrease  is  connected 
with  the  braking  of  CJ  head  part  because  of  magnetic  pressure  and  viscosity  forces’  work. 

The  work  carried  out  against  the  force  of  the  magnetic  pressure  causes  Lorentz  force  operation  with 

density: 

FL=jzBr=a(E2+Vr£,)Bt, 

where: 

Vr  radial  component  of  CJ  velocity  near  the  contact  surface; 

Ez  axial  electric  field  strength; 

<7  CJ  material  conductivity; 


pJ2 

2jir 


As  current  jz  induces  and  closes  in  skin-layer  limits  the  total  current  is  equals  zero  and 
Et~0  (accurately  Ez  determines  by  the  main  current  flow  conditions  and  this  value  can  be  neglect). 

Stationary  Navier-Stokes  equation  for  CJ  radial  velocity  Vr  on  arbitrary  radius  r  in 
cylindrical  coordinates  can  be  assumed  as  follows: 

(VV X=-— ^  +  «dKr, 

Pc  or 

where: 


where: 

P  pressure  in  CJ; 

pc,V  CJ  material  density  and  viscosity; 

R,r,r0  external,  current  and  initial  CJ  radiuses; 


A3.1 


A3 .2 


A3. 3 
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Vz  axial  jet  velocity. 

Axial  coordinate  z  =  0  corresponds  to  CJ  cross-section  near  the  internal  electrode,  z-l 
a  position  near  the  contact  surface. 

dV  dV 

In  assumption  that  — -  »  — —  Navier-Stokes  equation  transforms  to: 
dz  dr 

T,  8Vr  1  dP  d2Vr 

dz  p  dr  dz 2 

A3. 4 

Equations  A3. 1 ..  .A3. 4  can  be  solved  taking  into  account  the  following  assumption 

v,  =  K(r,z),  Vz  =  Vz(z),  P{z)=comt. 


Boundary  conditions: 


K\„  -o,vr\z__,=o 


V  -  V 

Y  z  z=o  0 


where: 


R  internal  CJ  radius; 

/  CJ  part  of  length  between  the  internal  electrode  and  the  contact  surface. 

Taking  the  integral  of  equation  (A3.4)  it  is  possible  to  receive  the  following  estimation: 

v  =_Wj_L+?±. 

dz  2  r 

Boundary  condition  Vr\r=0  =0  brings  the  relation: 

y  - 

dz  2 

A; 

In  equation  (A3. 5)  Vz  (z)  value  is  replaced  by  its  average  mean  Vz : 

tT  dVr  1  dP  d2Vr 

dz  p  dr  dz2 


Boundary  conditions: 
z  =  0,V~=o(^=o\ 


z  -l  ,Vr  =  0  for  viscous  liquid  on  solid  surface. 
The  solution  of  (A3. 6)  is  the  following: 


„  ldPl-z  Vz 

K  =~r - >  a  =  ”,  J]  =  pv 

tj  dr  a  v 
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From  equations  (A3. 6)  and  (A3 .7)  follows: 

dP  dVz  r  a 

dr  dz  2  l-z 


Equation  (A3. 8)  has  the  solution: 


n/x  dVz(R2-r2)  a  _ 

P(r)=-1-tr^  L~, —  +  PM. 

dz  4  l-z 


where: 


PM  magnetic  pressure  on  CJ  boundary  r  =  R  . 

Averaging  P(r)  over  the  radius  and  taking  into  account  thatP  -P0- 

( P{z )  =  const),  it  is  possible  to  receive: 

/  x  1  a  cfy  1  fjJ 2  Vz  Vz 

a°{,-y)=-6T,—2±-y(a°=J^x’y=r;'’  T) 


The  integral  of  (A3. 9)  with  boundary  condition  y{6)  =  1  is: 


y(z)=7r = exp{~  m%2\  m  =  L5{~ 
Vo  Vo) 


L)  Et lL  p-L 
J  Pcy  i 


Whence: 


y  =  y-  =  l  yfflf  =  J  exp{-  =  J 


For  *Jm  »  /  we  have  efr4tn  ~  1 .  Then: 


1  ( r  Y  P 
LL  G 

2{l)  PM 


A3. 10 


A3. 11 


A3. 12 


y  vz  U  Lm  \r0  J  Pc 


where: 

l  =  d  +  L 

m 

d  interelectrode  spacing; 

L0,Lm  penetration  deepness  without  and  with  current. 


A3. 13 


For  /  =  /  =  d  we  have: 

nun 
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( l  "1 

r 

7  d 

>  1  at  —  >  - 

KPm  J 

ro  1 

1  Pa 


2  R 


M 


A3. 14 


\  P 

At  J  =  400  kA,  r0~  1  mm,  — —  =  125 ,  d  >  10  mm. 


1  M 


In  case  L  »  d ,  l  ~  Lm  and  equation  (A3. 13)  has  a  view: 


f  T  ^ 


_ 


=  1.25 


\roJ 


C  p  ^ 

M 


\Pa  J 


(model  JVs  2), 


A3. 15 


For  CJ  with  r0  ~  1.5  mm,  L0  =  250  mm  at  current  J  ~  400  kA  factor - 5 .  At 

L 


current  J  ~  200  kA  the  factor  — —  ~  3  . 

Lm 

In  the  result  of  the  estimations  it  is  possible  to  note  the  following: 

•  The  model  of  CJ  electrodynamics  braking  is  developed  taking  into  the  viscosity.  However  the 
penetration  deepness  does  not  depend  on  viscosity. 

•  The  penetration  deepness  decrease  factor  is  determined  by  two  parameters  -  relation  between 
penetration  deepness  without  the  current  and  initial  CJ  radius  and  by  the  ratio  of  the  magnetic  pressure  to 
the  hydrodynamics  pressure. 

•  Calculated  penetration  deepness  decrease  factor  is  close  to  experimentally  observed  values. 
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___ _  APPENDIX  4 _ _ 

«PLANE»  MODEL  OF  JET  ELECTRODYNAMICS  BRAKING  INFLUENCING  BY  JET 

TRANSVERS  DEFORMATION 

The  hydrodynamics  model  for  which  the  assumption  of  pressure  constancy  along  CJ  is  not  needed 
(as  it  is  assumed  in  model  JV°  2)  is  discussed. 

The  mathematical  model  is  based  on  the  analytical  solution  when  the  viscosity  term  in  Navier-Stokes 
equation  is  neglected. 

Mote:  the  possibility  of  this  approach  was  looked  in  Appendix  3,  where  it  was  shown  that 
penetration  deepness  does  not  depend  on  CJ  viscosity. 

The  problem  can  pose  as  follows. 

It  is  necessary  to  determine  CJ  velocity  valuation  and  hydrodynamics  pressure  valuation  at 
axial  symmetrical  moving  at  the  conditions: 

•  At  r  =  0  CJ  is  asymptotically  represented  by  a  cylinder  of  r0  radius  the  axis  of  which  agrees  with  Z 
and  its  axial  velocity  is  V0 . 

•  At  r  =  / , infinity  target  is  pointed.  The  contact  surface  velocity  is  U  . 

•  In  CJ  plane  impact  there  is  formed  a  wave  of  radial  deformation. 

d  / 'JJ2 

•  Current  J0  flows  on  CJ  free  jet  surface.  This  current  generates  the  magnetic  pressure  PM  =  2~2-, 

where  R  -  distance  from  CJ  axis. 


Fig.  21.  Cylindrical  coordinate  system  connected  with  CJ  (Z-axis  is  perpendicular  to  the  plane 

of  the  drawing  in  0  point) 
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The  plane  option  of  the  problem  is  considered  in  spite  of  the  cylindrical  coordinates  are  selected  (see 
Fig.  21).  The  motion  is  considered  in  plane  with  the  coordinates  r,z,tp.  Axis  Z  is  perpendicular  to  the 

drawing  plane.  The  motion  along  Z  axis  is  absent,  but  Vr  =  Vr(r),  =  V,(r,<p). 

Two-component-stea.dy-sta.te  Navier-Stokes  equation  are  taken  into  account  [5]: 

(wK-^=--f, 

r  p  or 

(vvK+^=-4f 

r  P  dtp  A4i 


For  Vf  «  Vr  it  is  possible  to  receive: 


vdK=_l& 

dr  p  dr' 

W,  ,  l  SP 


\  dr  r  J  pr  dtp 


Continuity  equation: 


je(rK),isrr  o 

r  dr  r  dtp 


has  the  following  solution  (  Vp  «  Vr ) 

In  case  when  CJ  is  moving  as  a  whole 

M(tp)  =  where  y/  =  <pmax- 

Equations(A4. 1..A4.2),  give  the  result: 

r  dr  p  dtp 

From  the  equations(A4.3..A4.4)  neglecting  the  pressure  one  can  receive: 


The  integration  gives  the  following  relation: 


_d_f  Vi 

dtp\  2 


=  /W- 
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Having  Vr  -  Vr  ( r ),  one  can  find: 


K) 


=/W- 


Equations  (A4.6)  and  (A4.7)  are  identical  in  the  case: 

J_dP_ 
Pd(p' 


Integral  of  (A4.9)  is: 

f(y)v  =  -\—^~d(P  -  lPM~p(P)]=  [fa  )“  Po\  =  ~Pm  . 
J0  P  d<P  P  P  P 

where: 

hydrodynamics  pressure. 

Integral  of  (A4. 10)  if  to  take  into  account  (A4.9)  is: 

y  =  1  fMr  _  1  PMr 

9  2  M(y/)  2pM(y/)y/' 


The  radial  velocity  valuation  Vr  is  determined  from  equation: 


dP  _  dP_y  ~ 
dt  dr  r  9 


Or 


dy/ 

dr 


=  V„ 


From  (A4. 1 1)  and  (A4. 12)  taking  into  account  the  flow  conservation  equation  in  a  form 


A4.8 


A4.9 


A4.10 


A4.ll 


A4.12 


V  R  =  V  r 

v  rJX  y  0r  O' 


A4.13 


it  is  possible  to  receive 


l  dy  _lpM  l2 


1  +  r*—^-  =  -^-^Tr*,r*=j 


y/ dr*  2  PG  r0 


Vr  Vr  /  \ 

The  solution  of  (A4. 14)  at  condition  Vr  =  ■  =  0  0  and  Vr  [r0  )=V0  is: 


R  y nr' 


T  =  exP 


V  2  Pa  ro  J 


A4.14 


A4.15 
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Averaging  Vr  on  r*  in  interval  [o,/],  it  is  possible  to  receive: 

y  =  Z-  =  Ur  >)dr  *  =  f  ‘-JjJdf)  =  PL-Ler/J4  -  7E' ; 

^  17  ^  J  ^  ^  7  ,/^T  J  7  _fZ 7’ 


0  0 


2 


2  Vw 


(w  >  7,er/ Vw  *  ;),  4m  =  —  =  0. 7 

'  '  \2\PG  r0  \PG  r0 


Whence 


i  IPm  i 


v.  \  2  M  PG  r0 


A4.16 


A4.17 


For  /  =  lmin  -d  (I  =  d  +  Lm),  where  Lm  -  penetration  deepness,  the  condition  —  >  1  is  realized 

7. 


<7  , 

at  — >  J— 
r, 


W  \EL 


For  (J  =  ¥00  /c4,  r0  ~  1mm,  — 72J,  <7  >  15mm). 


M 


At  last  from  (A4. 17)  we  can  find:  ~  =  ^=r  =  —  =  —B- 

vr  U  Lm  nypa  r0 


or 


\ 

o 

a* 

1 _ 

4 

1 

_ 1 

Lm 

a? 
_ 1 

1 

1 _ 

(model  1) 


A4.18 


For  r0  ~  1mm,  L0  &  250  mm,  discharge  current  J  -  400  kA,  decrease  factor  - =  3. 3. Taking 


JUq 

for  the  discharge  current  a  value  J  =  200  kA  we  find  decrease  deepness  factor  to  be  equal  to  ——  =  2.3 . 
For  the  case  when  the  penetration  finishes  in  the  front  of  the  current  pulse  we  can  receive: 

"7(0  ^ max  .  *7 max  rp 

‘max  1 BH 

The  approximate  expression  for  penetration  deepness  is  followed  from  (A4. 18): 


J—  f  [E  )P¥. J—y  tdt  ( pmax  =p  (j  ))  Yil  =  L[E  K. 

PG  tbh  x  ’u  3ll*i  pa 


rM 

^ m 

\  Pa 

h 

—  ~0.6 
L 


pmax 

rM 

4 

IV 

L  pc  \ 

Jo  . 

A4.19 
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